Topological delocalization of two-dimensional massless Dirac fermions 
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The beta function of a two-dimensional massless Dirac Hamiltonian subject to a random scalar 
potential, which, e.g., underlies theoretical descriptions of graphene, is computed numerically. Al- 
though it belongs to, from a symmetry standpoint, the two-dimensional symplectic class, the beta 
function monotonically increases with decreasing conductance. We also provide an argument based 
on the spectral flows under twisting boundary conditions, which shows that none of states of the 
massless Dirac Hamiltonian can be localized. 

PACS numbers: 72.10.-d,73.21.-b,73.50.Fq 



The single parameter scaling theory of Anderson local- 
ization [31 predicts that the quantum transport of non- 
interacting disordered conductors is characterized by the 
beta function 



ding 
rflnL' 
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which encodes the variation of the dimensionless conduc- 
tance g with respect to the system size L. Once the 
value of the conductance at some length scale is known, 
the quantum transport at all length scales is constructed. 
The property of the beta function depends on the di- 
mensionality, and also on the symmetry class of the mi- 
croscopic Hamiltonian, such as spin rotation and time- 
reversal (TR) symmetries. Q In addition, the topological 
nature of wavefunctions also has a significant effect on 
quantum transport. 

In this paper, we discuss the problem of Anderson lo- 
calization for the two-dimensional (2D) two-component 
Dirac Hamiltonian subject to a random scalar potential. 



U = -ihvFcr ■ V + V{r). 
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Here, r G 



denote the standard Pauli matrices, 



and vp the constant velocity. The details of the random 
scalar potential V{r) will be specified later. 

The random Dirac Hamiltonian ([2]) is of direct rele- 
vance to the quantum transport of disordered graphene. 
[3| Although the band structure of clean graphene has 
two flavors (valleys) of two-component Dirac fermions, 
the intervalley scattering is rather weak since spatial pro- 
file of disorder in graphene is supposed to be smooth 
on an atomic scale. [5|, |6| A two-component single-flavor 
Dirac fermion can be realized, without doubling, on a 
surface of a three-dimensional Z2 topological insulator. 

The properties of the eigen functions for the ideal Dirac 
Hamiltonian (Eq. ([2]) without V) are well-known: The 
degeneracy point in the momentum space serves as a 



Dirac monopole for the Berry connection and wavefunc- 
tions in the momentum space pick up a tt p hase shift 
when transported around the Dirac cone. [lOj 

From the symmetry point of view, the random Hamil- 
tonian ([21) belongs to the symplectic symmetry class, as 
it possesses an "effective" TR symmetry 



(3) 



ll| The beta function of the 2D symplectic class shows 
the weak-antilocalization for large q, and there is a metal- 
insulator transition at g* ~ 1.4. [la. [isl. 

Although being a member of the symplectic symmetry 
class, there is growing evidence that the beta function of 
the random Dirac Hamiltonian ([2]) is qualitatively differ- 
ent from the conventional one for the 2D symplectic class: 
(i) Localization of non-rclativistic electrons for strong 
disorder can be understood in a picture in which bound 
states localized at potential minima overlap with each 
other. However, a Dirac fermion cannot be trapped by a 
potential well irrespective of the well depth [l5|, ll6[ , and 
hence is naively expected to have a strong tendency not 
to be insulating. This makes a physical picture for the 
strongly disordered regime of the Dirac fermions different 
from the conventional case, although physics of Anderson 
localization cannot fully be understood in terms of poten- 
tial trapping, (ii) As observed by Ando et ai, the Berry 
phase TT that is accumulated around the Dirac cone in the 
momentum space leads to a destructive interference be- 
tween a back scattering process and its TR counterpart, 
leading to the complete absence of back scattering. 
(iii) The non-linear sigma model (NLcrM, a field theory 
for diffusion modes) for the random Dirac Hamiltonian 
([2]) has a Z2 topological term. [13, 0, [3] It has a little 
effect in the metallic regime, but should change the renor- 
malization group flow in the strongly disordered regime. 
Ostrovsky et al. [l3| conjectured the NLcrM with the 
topological term has three flxed points (metallic fixed 
point, metal-semi-metal transition, and semi-metal at- 
tractive fixed point), (iv) There are numerical studies 
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FIG. 1: (Color online) The beta functions of the random 
Dirac Hamiltonian (closed circles) and the random spin- 
orbit model ([Sjl (open circles). The broken line represents 
the one-loop beta function of the conventional 2D symplectic 
class, f3{g) ~ l/rvg. [H 



that indicate the increase of the conductance with sys- 
tem size even for g < 1.4. 
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The purpose of this paper is to compute the beta func- 
tion of ^ numerically, and compare it with conventional 
system of the 2D symplectic class. We find that the beta 
function of the Dirac model is always larger than or equal 
to zero, showing that all states are delocalized even in 
the strong disorder regime, in contrast to the conven- 
tional case. We also provide a spectral-flow argument 
that clearly shows that the localization of Dirac fermions 
is forbidden. 

We compute the diagonal conductance (conductivity) 
of the random Dirac Hamiltonian ^ by evaluating the 
Kubo formula 
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where v is the velocity operator, v = i[T-l,r]/H = vpcr, 
f{E) is the Fermi-Dirac function at zero-temperature, 
r] in the energy denominator is a smearing factor, and 
|n) denotes an eigenstate with energy En of the Dirac 
equation in the presence of the random potential. 

The massless two-component Dirac equation cannot 
be regularized, without breaking the TR symmetry, by 
putting the system on a lattice. We thus work in the 
momentum space by introducing a hard cutoff at a suffi- 
ciently large momentum A. The eigcnstates {\n)} and 
energies {En} are then obtained by numerically diag- 
onalizing the Dirac Hamiltonian with disorder in the 
momentum-pseudospin basis. Typically, we take A ~ 
20 X 2tt/L, where about 2000 /c-points are included. We 
assume that the disorder potential is sufficiently weak so 
that the level broadening caused by disorder around the 
Dirac point is much smaller than the cut-off energy oc A. 

The smearing factor ij in the denominator of Eq. ^ 
accounts for the finite switch-on time of the electric field 
required for a dissipative current response. Physical ar- 



guments suggest that rj has to be at least as large as H/Tl 
where is the escape time from the system of interest. 
The escape time can be estimated from the Thouless 
energy (AE) by the uncertainty relation {AE)Tl — h, 
where AE is the eigenvalue difference between periodic 
and antiperiodic boundary conditions and_( ) is the geo- 
metric mean over disorder realizations. (23 . [2^ Indeed, 
it is reported in Ref. [20| that g is reasonably insensitive 
to 7] when rj ~ (AE). 

We assume that the scalar potential disorder V{r) is 
generated by randomly distributed impurities centered at 
R/, each of which contributes to V{r) with a scattering 
potential U{r — R/), 
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U{r - Ri) 
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potentials 
f/(q) ^ = 
potential. 



We considered two types of scattering 
U{r): the Gaussian correlated potential, 
w exp(— (7^Zq/2). and the Thomas-Fermi 
U{q) = u/{q + l'^^), where U{q) is the Fourier transform 
of C/(r), u represents the disorder strength, and Iq the 
range of the potential. Typically 5000 disorder configu- 
rations were used for averaging. The conductance was 
calculated for various sets of parameters, Ni, u, Iq and 
filling [Ep]', The number of scatterers iV^ was 1-10 times 
as large as the maximum number of carriers at each size; 
The range of the potential was changed upto 1 /30 of the 
minimal system size. 

We note that typical length scales are hardly de- 
termined from naive considerations at the Dirac point 
{Ep^ oo). Indeed the mean free path at the Dirac 
point, estimated by the golden rule, diverges for uncorre- 
latcd short-range scattering [Iq ^ 0) , while it vanishes 
for long-range Coulomb scattering (Iq — > oo)[2^. Never- 
theless we do not need the specific length scale since the 
beta function is defined as a logarithmic derivative in Eq. 

©• 

To compare our results with the conventional 2D sym- 
plectic class, we compute, by the same method, the beta 
function of the random spin-orbit (SO) coupling model 
given by 

(-iftV)V2m + y(r) 



H = 



-i{A(r),-*V} 



X (T 



(6) 



Note that the velocity operator in this model is spin- 
dependent. We assume an uncorrelated short-range dis- 
tribution for A(r) and V{r). 

Figure 1 shows the beta function of the Dirac model 
(filled circles) and the random SO couphng model (open 
circles). The latter agrees with the known behavior of 
the beta function of the 2D symplectic universality class: 
there is a metallic phase with weak anti-localization effect 
{(3{g) ^ I/ttp) when g is large, whereas there is a local- 
ized phase for small g; there is a metal-insulator transi- 
tion at g* 1.5 that separates the two phases. 
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(a) massless Dirac model 



(b) random SO model 





FIG. 2: (Color online) The evolution of energy spectra as a 
function of the twist angle (j> for the random massless Dirac 
model (left) and the random spin-orbit model (right). 



For large g, the beta function of the random Dirac 
model behaves similarly as that of the random SO model; 
this is as expected since when g is large, the Z2 topolog- 
ical effect in the NLcrM is small. 

We observe that the single parameter scaling holds rea- 
sonably well in both models as shown in Fig. 1. [23| In a 
sharp contrast to the conventional case, the beta function 
of the Dirac model monotonically increases with decreas- 
ing g well below g* ~ 1.4: A 2D massless Dirac fermion 
cannot be localized by a random scalar potential. (25| 
The numerical beta function of the Dirac model is well 
fitted by the one-loop beta function of the symplectic 
class even in the strongly disordered regime g 1. 

The absence of localization in the Dirac model can in- 
tuitively be understood by examining the spectral flow 
induced by twisting boundary conditions. Let us con- 
sider a finite and disordered system described by Eq. 
and impose the boundary conditions in both x and y 
directions, with phase factors exp(«02:) and exp(z0y), re- 
spectively. For simplicity we set (jiy = and discuss the 
energy levels as a function oi = 4>- The TR symmetry 
holds at (/) = 0, TT, where exp(i(f>) is real, leading to the 
Kramers degeneracy. We assume the cutoff to be infinity 
(the effects of the finite cutoff will be discussed later). 
Fig. [2Ia) shows an example of spectral fiow obtained for 
the 2D Dirac model with a specific disorder configuration. 
An essential observation is that Kramers pairs always 
change their partners as the energy spectrum evolves 
from = to tt; if the energy eigenvalues {En} are paired 
as • • • , {En,En+i), {En+2,En+3), ' ' ' at (/) = 0, then they 
are paired as •• • , {En-i, En), [En+i, E„+2), ■ ■ ■ at (/) = tt. 
Here eigenvalues En are ordered in ascending order [2^. 
In contrast, the non-relativistic electron system with SO 
coupling has a different type of the 'band-line topology' 
as shown in Fig. [2I^b); energy eigenvalues do not change 
their partners as the spectrum evolves from (f> — to tt. 

We can find the origin of this topological structure 
in the ideal spectrum. In the absence of disorder, the 
Dirac model has a set of eigenvalues En^.n ,s{4') = 



1/2 

{2ir/L)hvps [{ux + (j))^ + riy'j , where s = ±1 and show that (i) the conductivity hardly changes in a wide 



nxjUy 6 Z. For example, two degenerate states at = 
with zero energy (s = ±1 and Ux = riy = 0) become apart 
as (j) increases and never stick together; each couples with 
other partners at = tt. As we introduce disorder, en- 
ergy eigenvalues move around but the way eigenvalues 
are paired between = and tt can never be altered, 
since each Kramers doublet remains sticked at = and 
TT. In other words, it is impossible to change the topology 
of the 'band-line' continuously from the type of Fig. ^a) 
to (b), without breaking the TR symmetry. 

If a state is exponentially localized, its eigen energy 
must be insensitive to the boundary phase factor, i.e., 
the 'band width' of En{4>) is exponentially small com- 
pared with the average level spacing In the Dirac 
model, however, it is impossible because all the band lines 
are connected through the Kramers doublets at = 0, tt 
so that the band width cannot be smaller than the level 
spacing. We thus conclude that there are no localized 
states in the Hamiltonian In the non-relativistic elec- 
tron system, in contrast, the structure of the spectrum 
in Fig. [21(b) does not, at least, prohibit localization, and 
states indeed tend to be localized for strong disorder. 

We note that, in order for the above argument to be 
valid, we have to assume that the energy band continues 
from —00 to 00. Indeed, if we have a finite cutoff, the TR 
symmetry must be broken either at = or tt. Although 
this may alter the band-line topology around the band 
edges, the low-energy states around the Dirac point are 
hardly affected as long as the disorder potential is long- 
ranged and the cut-off is large enough. With increasing 
the disorder strength, one would naively expect Anderson 
localization first takes place at band edges (cut-off) and 
the Dirac point (can be viewed as a point at which two 
band edges meet accidentally) . The former goes away as 
we send the cutoff to infinity, while the latter is protected 
from localization by the topology of the spectral flow. 

Although the honeycomb lattice system involves a cou- 
pling of the two valleys (flavors), a similar delocalization 
effect should manifest itself when intervalley scattering is 
negligibly weak. On the other hand, when atomic-scale 
scatterers dominate, the intervalley scattering random- 
izes the Berry phase and the nature of interference is 
changed to enhance localization. In the Dirac band, 
the inter-valley scattering time depends on the Fermi en- 
ergy as oc l/|i?F|01j and thus is more important in the 
highly doped regime. 

The present calculation suggests that the Dirac 
fermion system exhibits the positive magnetoresistance. 
On the other hand, the recent graphene experiments in- 
dicate somewhat complicated situations: A magnetore- 
sistance study [i^l clarified that highly doped epitaxial 
graphene exhibits a crossover between positive and neg- 
ative magnetoresistance induced by changing the tem- 
perature as expected theoretically in Refs. p, [28[. For 
isolated single graphene sheetsQ, however, experiments 
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range of temperature near the Dirac point, while (ii) 
the magnetoresistance is weakly positive at low carrier 
densities. [2§| Although a number of theoretical scenar- 
ios have be en p roposed, including effects of microscopic 
ripp les [29I . [30| . trigonal warping terms [2§|, and edges 
[3l|, there is no consensus at this moment. Taking into 
account these effects in addition to the random scalar 
potential in Eq. ^ will be done elsewhere. 

Although our focus in this paper is on 2D, the argu- 
ment based on the topology of the spectral flow applies 
equally well to the 1 and 3D two-component massless 
Dirac fermion with the effective TR symmetry: A two- 
component massless Dirac fermion cannot be localized by 
a random scalar potential in all 1, 2, and SD.js^l Since d- 
dimensional two-component massless Dirac fermion can 
be viewed as a gapless boundary mode of Z2 topological 
insulators in {d+l)D {d = 1, 2, 3)0, H, 0], our discussion 
above concludes that a surface of a (strong) Z2 topolog- 
ical insulator is always metallic, robust against disorder. 
This is consistent with the speculation in Ref. 0] in the 
context of quantum spin Hall effect. Such metallic sur- 
face states can be called a "topological metal" . 

After completion of this work, we became aware of 
a similar numerical result at the Dirac point, obtained 
independently in Ref. 
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